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If the expression of the topological charge density operator, suggested by fermions obeying the Ginsparg- Wilson 
relation, is employed, it is possible to prove on the lattice the validity of the Witten-Veneziano formula for the 
rj' mass. Recent numerical results from simulations with overlap fermions in 2 (abelian Schwinger model) and 
4 (QCD) dimensions give values for the mass of the lightest pseudo-scalar flavour-singlet state that agree with 
theoretical expectations and/or experimental data. 



1. Introduction 

In a nut-shell the C/a(1) problem lies in the fact 
that, in the absence of ?7yi(l) anomaly contribu- 
tions, the unphysical bound r7ioctot,7=o < V^mT^ 
holds §. In 1976 't Hooft §1 pointed out that 
the resolution of this problem had to be related 
to the existence of topologically non-trivial gauge 
field configurations in Euclidean QCD. Dwelling 
on this observation, it is possible to obtain an ex- 
plicit formula for the ry' mass either in the 't Hooft 
limit {Nc^oo, with g^N^ and Nf held fixed g), 
as was done in [Q, or by assuming that anoma- 
lous flavour-singlet axial WTIs retain their valid- 
ity order by order in an expansion in u = Nf /N^ 
around m = 0, as argued in In both cases one 
can derive the "leading-order" Witten-Veneziano 
(WV) relation 

= I d^x{Q{x)Qm , (1) 

' J YM 

where F^r ~ 94MeV is the pion decay con- 
stant, A the "topological susceptibility" and Q = 

2 

^:^etj.iypatT[Ff,^Fpcr] the topological charge den- 
sity. The notation (...)|ym means that the QQ- 
correlator is to be computed in the pure Yang- 
Mills (YM) theory, i.e. in the absence of quarks. 

The formal relations in eq. (|l|), when translated 
in any regularized version of QCD, such as lattice 
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QCD, become more involved and quite a number 
of subtleties have to be dealt with in order to be 
able to determine the correct field theoretical ex- 
pression of A, which should be used in the above 
equations. Two problems need to be solved to 
make the formulae in (]l|) precise and of practical 
use. One has to i) find a properly normalized lat- 
tice definition of the topological charge density, 
Ql] ii) subtract from the Ql{x)Ql{0) product 
appropriate contact terms, as required to make it 
an integrable (operator-valued) distribution. 

A rigorous derivation of a formula for the rj' 
mass, which could be unambiguously employed 
in numerical simulations, can be given if use 
is made of the (lattice) regularized anomalous 
flavour-singlet axial WTIs of full QCD in order to 
construct a properly defined (QlQl) correlation 
function. The remarkable result of this analysis 
is that no subtraction is needed to this end in the 
chiral limit, if the definition of lattice topological 
charge density suggested by fermions obeying the 
Ginsparg-Wilson (GW) relation (e.g. overlap 
fermions [^) is employed. 

2. GW fermions and the rj' mass formula 

Regularizing the fermionic part of the QCD 
action using GW fermions offers the great 
advantage over the standard Wilson discretiza- 
tion ^ that global chiral transformations can 
be defined, which are an exact symmetry of the 
massless theory, as in the formal continuum limit. 
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This follows from the relation 75!) + D75 = 
D^bD 0, obeyed by the GW Dirac operator, 

In this regularization the Ua{^) anomaly 
is recovered a la Fujikawa [0, because of the 
non-trivial Jacobian that accompanies the change 
of fermionic integration variables induced by a 

lattice transformation. As a consequence 
the (anomalous) flavour-singlet WTIs in the pres- 
ence of Nf massless fermions take the form 



2N 



where O is any finite (multi)local operator, 6^0 is 
its local variation and V^A°(x) is the divergence 
of the singlet axial current. 

Neither Tr[j^D{x,x)] nor A^{x) are finite op- 
erators, but finite linear combinations, Q and Aj^, 
can be easily constructed pH by writing 

Q{x) = ^Trh5D{x,x)]-^V^Alix) (2) 

i°(x) = (l-Z)4(x). (3) 

where Z is logarithmically divergent to lowest or- 
der in perturbation theory and vanishes as u ^ 0. 

With the above definitions the renormalized 
singlet axial WTI becomes 



V,{A»0) = 2Nj{Qix)d) + {510) 



(4) 



This equation shows that Q and A° are correctly 
normalized and that there is no power-divergent 
mixing of Tr[75_D(x, x)] with the pseudo-scalar 
quark density, which would bring in a dangerous 
liner aly divergent mixing coefficient. 

A formula for the 77' mass can be obtained from 
the above WTIs, by observing that the rj' mass 
must vanish as u for massless QCD to have 
no ^-dependence. For that one starts by defining 
the lattice Green function |^ 

Xtdp) ^T^r f d^xe-^P^W,{Al{x)Qm , (5) 
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^We set the lattice spacing, a, equal to 1. In all lattice 
formulae the continuum limit is, however, understood. 
*A contact term, CT(p), should be added to the r.h.s. of 
eq. (^) to make it finite at p ^ 0. CT(p) is a polynomial 
of degree 4 in p, which vanishes at p = 0, because the 
r.h.s. of eq. is certainly finite (actually zero) at p = 0. 
CT{p) plays no role in the argument below, since we will 
be finally only interested in the value of XtL at p = 0. For 
brevity we will not indicate it explicitly in the following. 



with Q given by eq. (g) . In the full theory, where 
the r]' is massive, one has Xti(O) = 0. On the 
other hand, since, as we obseved above, at the 
chiral point, the ry' mass vanishes as m ^ 0, only 
the 77' pole will contribute to Xthip) in this limit, 
leading to the relation 



lim lim XtL (p) 

p— i-O u—*0 



2Nf 



u=0 



(6) 



Thel.h.s. ofeq. (|) can be evaluated by using the 
WTI (|) with O = Q. Recalling that Z vanishes 
when u ^ and 6^Q = 0, we obtain 



p2 
2Nf " 



ti=0 



lim lim 



;iTr[75i?(x,x)]iTr[75i^(0,0)]). 



(7) 



The limits u —^ and p in this expression 
can be readily performed in the order indicated, if 
one can assume that taking the first limit simply 
amounts to setting the fermion determinant equal 
to unity. One gets in this way 
i?2 



2Nr^' 



u=0 



d^x{-TT[j5D{x,x)]-TT[j,D{0,0)]) (8) 

2 2 YM 

The restriction to pure YM theory, indicated in 
eq. (^), is an obvious consequence of the fact that, 
for a Green function of only gluonic operators, ne- 
glecting the fermion determinant is tantamount 
to limiting the functional integral to the pure 
gauge sector of the theory. 

A relevant question at this point is to ask at 
which value of Nc eq. (||) is supposed to be valid. 
The answer depends on the behaviour of QCD 
with Nf. Various scenarios are envisegeable. 

1) In the most favorable situation, in which the 
limit u — !■ of Xthip) exists and is equal to the 
value it takes at Nf = (fermion determinant 
equal to 1), formula (^) is valid for any value of 
Nc and for each N^ it yields the mass of the 77' 
meson (at leading order in u = 0) in the world 
with the corresponding number of colours |^] . 

2) If quenching can be attained only in the limit 
in which the number of colours goes to infinity, 
then eq. (^) will yield a formula for the m^, valid 
up to 0[1/ Nc) corrections 0]. 
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3) Finally it may happen that taking the limit 
M — > does not correspond to quenching. In this 
case one cannot pass from the fairly complicated 
eq. (0) to the more useful formula (||). 

Remembering that ^Tr[75D(a;, a;)] should be 
identified with the topological charge density , 
we conclude that eq. (^) can be rewritten in the 
very suggestive form 



— ^1 

2Nf 



lim 



M = 



V 



(9) 



where {{nj^ — ni^)'^) is the expectation value of the 
square of the index of the GW fermion operator, 
D, and V is the physical volume of the lattice. 

In the form the r/' mass formula can be di- 
rectly compared with overlap fermion simulation 
data. Existing numerical results in 2 (abelian 
Schwinger model [||) and in 4 (QCD di- 
mensions agree quite well with theoretical expec- 
tations and/or experimental numbers 

3. Conclusions 

We have shown that a formula for the rj' mass 
can be rigorously derived in lattice QCD, exploit- 
ing the anomalous flavour-singlet axial WTIs of 
the theory. If fermions obeying the GW relation 
are used, there exists a natural definition of topo- 
logical density which is correctly normalized and 
for which the naive form of the WV formula holds 
without the need of introducing any subtraction. 

In the literature two sorts of approaches have 
been proposed to deal with the problem of com- 
puting A on the lattice, which have led to nu- 
merical values as good as those obtained with the 
present method Q . The first one is based on a di- 
rect field theoretical definition of A ||T^ that takes 
into account the need for the renormalization of 
Ql and the subtraction of the operators and 1 
in the short distance expansion of QlQl- The 
second one makes use of the notion of "cool- 
ing" ||l5| to carry out the necessary operations of 
renormalization and subtraction. Both methods 
are reliable to the extent that they are able to cap- 
ture the topological properties of the gauge field 
configurations that determine the number of zero 
modes of the Dirac operator. Simulations based 
on the geometrical definition of of rcf. [p6[ 
have not yet led to comparably good results [|17|. 
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